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Abstract 

Communications in dispersive direct-sequence code-division multiple-access channels suffer from intersymbol 

, and multiple-access interference, which can significantly impair performance. Joint maximum a posteriori probability 

equalization and multiuser detection with error control decoding can be used to mitigate this interference and to 

^ , achieve the optimal bit error rate. Unfortunately, such optimal detection typically requires prohibitive computational 
OO . 

' complexity. This problem is addressed in this paper through the development of a reduced state trellis search detection 



o 
in 
o 



X 



algorithm, based on decision feedback from channel decoders. The performance of this algorithm is analyzed in the 
large-system limit. This analysis and simulations show that this reduced complexity algorithm can exhibit near-optimal 



C/2 ■ performance under moderate signal-to-noise ratio and attains larger system load capacity than parallel interference 
cancellation. 



I. Introduction 

Over the past two decades there has been considerable research on direct sequence code division multiple 
access (DS-CDMA) communications, which offers the advantages of soft-capacity limit, inherent frequency 
diversity and high data rate ifTSl . and which is the fundamental signaling technique of third-generation 
(3G) cellular communications and other emerging applications.. However, DS-CDMA suffers from inter- 
ference, particularly multiple-access interference (MAI), which is due to the non-orthogonality of different 
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users' spreading codes, and intersymbol interference (ISI), which is caused by multipath fading in high-rate 
systems. It is well known that mitigation of these types of interference can substantially improve system per- 
formance. In recent years, much progress has been achieved toward such mitigation through the application 
of equalization (EQ) and multiuser detection (MUD). 

Equalization can be roughly categorized into three classes. One class is based on maximum likelihood 
(ML) detection ifTSl . which can be implemented efficiently with the Viterbi algorithm (VA). A second 
class is linear equalization based on some criterion, such as minimum peak distortion or minimum mean 
square error (MMSE). The third class is decision feedback equalization (DEE), in which previously detected 
symbols are used to cancel the intersymbol interference lITSl . In multiuser detection, we can find the 
counterparts for these three kinds of equalizers. In particular, ML based multiuser detection and MMSE 
multiuser detection are both well known [|24| . As a combination of both techniques, the problem of joint 
ISI an MAI mitigation is discussed in 

In recent years, the turbo principle, namely the iterative exchange of soft information among different 
blocks in a communication system to improve the system performance, has been applied to equalization 
and multiuser detection in channel coded systems, thus resulting in turbo equalization |JTll and turbo mul- 
tiuser detection ll23l . In these algorithms, soft decisions from channel decoding are fed back to be used 
a priori probabilities by a maximum a posteriori probability (MAP) based equalizer or multiuser detector 
and enhance the performance iteratively. However, the computational cost of MAP based detection is pro- 
hibitive for large numbers of users or long delay spread. Thus, it is necessary to reduce the complexity 
of such iterative algorithms for practical applications. Eor memoryless synchronous multiaccess channels, 
the MAP turbo multiuser detection can be simplified to parallel interference cancellation (PIC) f\\, whose 
performance can be enhanced with an MMSE filter lE^ or a decorrelating interference canceller Ifim . For 
systems with memory, an alternative way to simplify these detectors is to reduce the number of states by 
either truncating the channel memory to a fixed order or eliminating the states with reliable decisions at the 
channel decoder. The former strategy has been used in the equalization of ISI channels |l6| Q lfT3l and 
in asynchronous multiuser detection of channel coded CDMA systems lfT6ll . while the latter scheme is of 
dynamic complexity and is widely used in iterative decoding algorithms. In speech recognition ll25l . this 
scheme is applied to trim the acoustic or grammar nodes with low metrics. In early detection based decoding 
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of parallel turbo codes the trellis is simplified by splicing the state with reliable a priori probabilities. 
Similar strategies have been applied to the iterative decoding of LDPC codes |,8J and concatenated codes 

lH. 

In this paper we study channel coded DS-CDMA systems operating over frequency selective fading chan- 
nels. Although joint detection and decoding can achieve higher channel capacity, we consider only systems 
in which detection and decoding are separate due to the increased feasibility of such systems in practical 
applications. In particular, we consider joint equalization and multiuser detection based on MAP detection 
with decision feedback (MAP EQ-MUD) from the decoder. Similarly to MAP turbo equalization or mul- 
tiuser detection, the MAP EQ-MUD also suffers from prohibitive computational cost. In this paper, we first 
decompose the multiuser trellis into single-user trellises, thus linearizing the number of states in terms of 
the number of users. Then we apply both the state-reducing techniques in channels with memory, namely 
shortening the channel memory to a fixed order and partitioning the decision feedback of channel decoders 
into unreliable and reliable sets of symbols using a simple confidence metric. The states are constructed 
with the unreliable set and the interference from the reliable set is cancelled using the soft decisions from 
the decoder. We call this reduced complexity algorithm reduced state equalization and multiuser detection 
(RS EQ-MUD). 

This paper is organized as follows. Section II introduces the system model, in which the signal model and 
decoder are explained. The optimal MAP EQ-MUD algorithm is described in Section III, while Section IV 
is focused on developing the RS EQ-MUD algorithm. An asymptotic analysis of the system performance is 
carried out in Section V, and corresponding numerical results are given in Section VI. Final conclusions are 
drawn in Section VII. 

II. System Model 

A. Signal model 

For a channel-coded DS-CDMA system, let K denote the number of active users, the spreading gain 
and M the coded symbol block length. Denote the symbol period and chip period by Tg and Tc respectively, 
and note that Tg = NT^. For user k, the binary phase-shift keying (BPSK) modulated continuous-time signal 



at the transmitter is given by 

M 

hit) = J2^k{i)skit - iT,) 

i=l 

M N 



n=l 



where bk{i) is the i-th (binary) channel coded symbol sent by user k, ^(t) is the chip waveform, and s^^\n 



The superscript i in s^*^ (n) 



is the normalized binary spreading code of user k, which satisfies 
implies that the spreading code varies with the symbol period, since long (aperiodic) codes are considered 
in this paper. 

The signal (1) passes through a frequency selective fading channel whose impulse response is hk{t), and 
the channel output is the convolution of the input signal and the channel response: 

rk{t) = fk{t)-khk{t); 

M N 

= J2^''('^J2^k\^)9kit-iTs-n%), (2) 

i=l n=l 

where 

gk{t - iT, - nT,) = ^{t - iT, - n%) * hk{t) . (3) 
Assuming synchronous transmission among users \ the signal at the receiver is given by 

K 

k=i 

KM N 

= J2J2^k{i}J2^k\^)9k{t-tTs-nn). (4) 

fc=l i=l n=l 

The received signal is sampled at the chip rate ^, and the corresponding discrete output of the sampler at 
chip period / is 

KM N 

k=l i=l n=l 
K M 

= J2J2bkii}ht\l-iN), (5) 

k=l i=l 

^Note that this assumption of synchronous transmission does not limit the generality of this model, since delay offsets between users can be 
incorporated into the channel impulse response hi{l), hK{l)- 



where 

h^\l)^s^\l)i.g,{lT,). (6) 

Suppose the support of h^^\l) is (0, (L — 1)A^), where L is the dispersion length. We can simplify (5) to a 
vector form. (For notational convenience, we henceforth use t to designate discrete time advancing at the 
symbol rate.) This vector is given by ^ 

y{t) = {y{tN),y{tN + 1), + 1)N - l)f , 



and 



(i) 



{j)^{h^\jN),h^\jN+l),...,h^\{j + l)N-l))y 



Using the vector form and considering the thermal noise at the receiver, we write 

K t 

yW = E E bk{i)h^\t-t)+n{t), t=l,2,...,M (7) 

k=l i=t-L+l 

where n(t) is additive white Gaussian noise (AWGN) that satisfies E {n(t)n(t)^ j = a'^I^^^. The devel- 
opment in the remaining of this paper will be based on this discrete model (7). 

B. Equivalent spreading code 

We term h^*^ [t — i) the equivalent spreading code of the i-th symbol of user k in the t-th [i < t) symbol pe- 
riod. In order to explore the properties of the equivalent spreading code, we need to place some assumptions 
on the discrete channel response gk{lTc) in (6): 

• Causality. gk{lTc) = when I < 0. 

• Normality. We assume that {gk{lTc)}i is a complex-valued Gaussian random sequence with zero mean 
and exponentially decaying variance 

4(1) = (8) 

where is the decay factor and satisfies e~'^'=^ pa and lim^v^oo Y^^i ^li^) = 1- 

• Independence. We assume that the scattering caused by fading is uncorrelated, which means 

E {gk{lTc)gn{jTcr} = 0, if j or k^n. (9) 
denotes transposition, H denotes conjugate transposition and * denotes conjugation. 
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Thus gk{lTc) and gn{jTc) are independent since they are jointly Gaussian. 
In non-dispersive multiple-access channels, the cross correlation of the spreading codes plays a key role 
in the system performance. Thus we need to discuss the cross correlation between the equivalent spreading 
codes. Define the equivalent cross correlation between the equivalent spreading codes of user k and user n 
to be 

prn(^,j) = (hi*^(0)''hM(j)- (10) 
Note that ^^'^(i, i) is a random variable since h|(i) and h^(i) are random. 

C. Channel decoder 

The diagram of the system discussed in this paper is given in Fig. 1. At the transmitter, the information 
symbols are encoded, interleaved with an infinite length interleaver and spread. The BCJR algorithm ||2|| is 
used in the channel decoder, which follows a deinterleaver, to obtain the a posteriori probability P{bk{t) = 
b\yf^), where yf = {y{t)\t = 1, M}. This probability can be decomposed into two parts, 

pm) = b\yr)^p{h{t) = b\ymi{b), 

where Cli^) is the extrinsic information about bk{t) from the other coded symbols. We use the extrinsic 
information to construct the soft decision feedback: 

hit) = 2(1(1) -1, (11) 

which is used to cancel both the ISI and MAI. The estimation error is denoted by Abfc(t) = bk{t) — bk(t); its 
expectation is zero due to symmetry and its variance can be obtained by simulation. 

III. Optimal MAP Turbo EQ-MUD 

In this section we develop the BCJR algorithm fT\ for MAP EQ-MUD in a way similar to that used for 
turbo multiuser detection in ll23ll . 

First we define the state at symbol period t as the set of all symbols of all users from symbol period 
t- L + 2tot: 

St = {bkim = t-L + 2,...,t, k = l,...,K}. (12) 



Thus for each symbol period we have 2^^~^^^ states with which to construct a trellis. We call state m and 
m' compatible when the state can transit from m to m' and denote this condition by m =^ m'. 
Next we define some intermediate variables 0: 

• forward probability: Q!t(in) = P(St = in,y*); 

• backward probability: Ptim) = P{yf^^\St = m); 

• transition probability: 7i(m',m) = P{y{t),St = m|St_i = m'). 

The forward and backward probabilities can be computed recursively via the equations 

at{m) = ^at_i(m')7t(m',m), (13) 
m' 

and 

A(m) = ^A+i(m')7<+i(m,m'). (14) 
m' 

Invoking Bayes' formula, we can rewrite the transition probability as 7t(m, m') = F(y(t)|Si = m, S^^i = 
m')P(St = ni|S(_i = m'). Since we assume that the interleaver has infinite length, the symbols of different 
users and different symbol periods are independent. Thus we have 

P{St = m\St-i = m') = P{h{t) = h,...,bK{t) = bK\bi,..jK em) 

K 
k=l 

= A,(m), (15) 

which can be regarded as the a priori probabilities of the symbols at symbol period t. When no a priori 
information is available, we assume the symbols to be uniformly distributed. Since the channel decoder can 
provide the extrinsic information from the other coded symbols, we can feed the soft outputs of the decoder 
back as the a priori probability. Hence 



if m' m. 



where yni,m'(^) is the estimated received signal due to the states m and m': 

K t 

ym,m'W = I] Yl ~bk(i}bk\t-i), em'Um. (16) 

fe=l i=t-L+l 

When the forward and backward probabilities are available, we can compute the joint probability 

P {St = m,yf ) = at(m)A(m), (17) 
and then compute the a posteriori probability: 

P (hit) - 6fe|yf ) « J] P (St = m, yf ) , (18) 
where oc denotes proportionality. The appropriate normalization is 

P {b,{t) = l\y^) + P {hit) = -l\y^) = 1. 

To avoid the reuse of the extrinsic information, we need to cancel the a priori probability. Therefore the soft 
input to the decoder, denoted by Ck{t), for symbol t of user k is normalized by the soft decision feedback: 



Ck{t) 



' p(6fcW=i|yf ) " 


/(bfc(t)=-i|yf)_ 




' C(fefeW=i) " 




Lc(&/c(t)=-i)J 



(19) 



IV. Reduced State EQ-MUD 
A. Independence assumption and trellis decomposition 

Facing the same problem as either equalization or multiuser detection, the optimal MAP EQ-MUD suffers 
from prohibitive computational complexity because the number of states increases exponentially with the 
number of users K and the dispersion length L. Thus, the optimal MAP EQ-MUD is primarily of theoretical 
value and cannot be implemented for many practical applications. 

In this paper we mitigate the complexity of the MAP EQ-MUD, by decomposing its trellis into single-user 
sub-trellises, thus linearizing the number of states with respect to the number of users. For user k, we define 
the sub-state at symbol period t as 



S^^^{hu{t-L + 2),...Mt)}- 
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The sub-states belonging to one user construct a single-user sub-trellis. Then any state defined in (12) is the 
union of the corresponding sub-states of all users. In order to distinguish sub-state and state, we use non- 
bold fonts to designate the sub-states in the rest of this paper. Similar to Section III, we define the forward 
and backward probabilities for the sub-state of user k: 

af(m) = P{S!^^m,y{), (20) 

and 

P'^im) = P {S'^ ^ m,y^^,) . (21) 

Compared with the analogous definition of Section HI, the definition of the backward probability Pt{m) is 
slightly different here; this will facilitate the application of Assumption IV. 1, which follows immediately, on 
the backward probability. If the BCJR algorithm can be confined to each sub-trellis, the number of sub-states 
will be reduced to K2^~^. However, the forward and backward probabilities defined in Section EI involve 
joint distributions of different users' symbols; thus the sub-trellis searches for different users are coupled, 
which prohibits the exact decomposition into single-user trellises. 

However, we can approximate the joint distribution with the product of marginal distributions, which 
results in the following assumption. 

Assumption IV. 1: For any state defined in (12) which is the union of the correspongding sub-states, 
St^ULiS^, we have 

K 

p^yi) = n^('^'iyi). (22) 

k=l 

and 

^(S*|y*^) = n^(^*lym)- (23) 

fc=i 

This assumption is based on the fact that P {S^\ y{ ) is concentrated around and 1 , provided that the noise 
power is small enough. It can be validated by the simulation results in Section VI, which state that, with 
Assumption IV. 1, the reduced complexity algorithm in this paper can achieve near optimal performance 
in moderate energy region. Under this assumption, which is only an approximation, we can decompose 
the forward and backward probabilities of a state into the product of the probabilities of the corresponding 
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sub-states: 

/O N _ nf=l O^KSt) .r,.. 

(^(yi)) 

and 

A(S.) ^ M) . 

We can neglect the common factors (^(yi))^^ and (F(y^i)) ^ ^, which do not affect the final result. 

With this assumption, we can develop recursive formulas similar to (13) and (14) with respect to the 
sub-states. In particular, for the forward probability, we have 

= E E ^(S. = m",S,_i = m',y*) 
m' mem" 

^ ^(yl|Si = m",Si_i=m')P(Si = m"|Si_i=m')P(Si_i = m') 

m' mem" 

- P(yr\S,_i=m')P(y(t)|S, = m",Si_i=m')P(S, = m"|Si_i = m') 

m' mem" 

= Ell n <iW)A,^(m',m)), (26) 
m' \ \i=i,m<em' / / 

where (m', m) is the probability of transition from state m' to sub-state m. 

mem" 



mem" 



\^ A (^"\_L_ ( lly(^)-^m",mK^)ir \ 

^ mS" ^Vi^^^' V ) ' '''' 



where Win",m'(^) is the estimated received signal due to the states m' and m": 

K t 



Wm",m'W = E E bn{i)y^fit-i), (28) 

n=l i=t-L+\ 



where 6„(i) e m' U m". 
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For the backward probability Pt{m), we can obtain a similar recursive formula with the same manipula- 
tions as the forward probability: 

Ptim) = E^(y*+2,St+i = in')P(y(t + l)|S, = m",S,+i=m')P(Si = m"|S,+i=m') 

mem " m' 

= Ell n /5mK))0t'+iKm')), (29) 
m' \ \i=i,miem' / / 

where ^^^.^(m, m') is the probability of transition from sub-state m to state m'. With the same manipulation 

as (m', m), we have 

mem'' V^^< V 2(7„ J 

where Vni",m'(^ + 1) is the estimated received signal due to the states m' and m": 

K t+l 

Vm'Mn'(i+l) = E E ^n(i)h«(^ + 1 - (31) 

n=l i=t-L+2 

where e m' U m". 

When the forward and backward probabilities are available, we can compute the joint probability in a way 
similar to (17): 

P(g>=„..yr) = y''1 - (32) 

P{S^ = m) 

And the likelihood ratio for input to the decoder can be computed with (18) and (19). 

B. Reduced state trellis 

With the above independence assumption, we have reduced the number of sub-states to K2^~^. However, 
the BCJR algorithm of each sub-trellis is still coupled to the other sub-trellises because of the existence of 
MAI. In (26) and (29) the number of terms in the double summations increases exponentially with K. Thus 
the computational complexity of the decomposed trellis still remains prohibitive for practical implementa- 
tions, and we need reduce further the number of sub-states. We can accomplish this with the aid of the 
decision feedback from the decoders. The underlying philosophy is to partition the decision feedbacks into 
reliable and unreliable sets, denoted by R and R. The trellis is constructed with the symbols in the unreliable 
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set while the symbols in the reliable set are considered to be known and are directly cancelled from the 
original signal. 

Suppose that the decision errors of the decoders are {A6fc(t)} ^ ^ ^ ^ ■ We sort the absolute values of 

i = 1,2, ...,M 

these errors and select the largest K ones to form the unreliable set R while the remaining symbols comprise 
the reliable set R. We define k = as the normalized search width, which is also the probability that a 
decision feedback symbol is unreliable. 

When we are constructing the sub-states at symbol period t, we consider only the symbols from t — Li + 1 
to t where Li is an integer with I < Li < L, thus shortening the channel memory in a way similar to 
||6| lfT3l . Li can be selected so that the decision feedback errors of the symbols earlier than Li compose 
only a small portion of the interference due to the exponential decay of the channel. With the unreliable set 
R and notation R(t) = {(A;, z) G R, t — Li + 1 < z < t, 1 < A; < K}, we define the sub-states of user k at 
symbol period t: 

St = {bumk,^) eK{t) ori = t]. 

Observe that after reducing the sub-states, the numbers of sub-states of different users may be different and 
depend on the selection of the unreliable set. Thus the total number of sub-states at symbol period t is 
Efc=i 2''*+\ where Hk is the number of unreliable symbols of user k from symbol period t — Li to t — 1 
and has the constraint Yl!k=i f^k = the size of set R(t). Assuming that fik is Bernoulli distributed, the mean 
value of the number of sub-states at symbol period tisK J2m=o (lI^i) '^'"(1 ^ «;)^^~^~'"2'"+^. 

The BCJR algorithm for decoding the reduced sub-state trellis is also confined to the unreliable set. The 
transition probability is nonzero only when m' ^ m and m' C R(t) for (26) and when m ^ m' and 
m' C R(t + 1) for (29). The recursive formulas are rewritten as follows: 

«tM = E (( n A,^m',m)j , (33) 

m'cR(t) V \^=l,"^,em' / / 

and 

A'M = (in AViK) 1 0tiK«i') I • (34) 
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And we rewrite (28) and (31) as follows: 



Wm'Mn'W= J2 bn{^H\t-^)+ ^n«h«(t-2), (35) 

and 



vm",m'(t + l)= k{^M\t + l-^)+ 6„(2)h«(t + 1 - (36) 

where &„(i) is the reliable decision feedback about symbol 6„(z) from the decoder. 

Figure 2 illustrates the search of the sub-trellis with K = 2 and Li = 2. The state transitions from 
symbol period t — 1 to t are labelled with arrows in the figure. Since 1)2 (t — 1) is reliable, the corresponding 
interference is cancelled directly. Thus we consider only the transitions from S^_-^ to and 5"^^ since 6i (t— 1) 
is unreliable. 

Since the decision feedback errors are unknown to the receiver, we cannot use {Abfc(t)} to construct the 
unreliable set and it can only provide an upper bound for the performance. It is easy to show that the larger 
the absolute value of bk{t) is, the more reliable the soft decision feedback is, when only the soft decision 
feedbacks are available. Thus we call |6| the confidence metric for symbol h and, use it to construct the 
unreliable set. 

C. Distribution of the unreliable set 

For asymptotic analysis of system performance in Section V, we need to obtain the asymptotic distribution 
of decision feedback error in the unreliable set. Suppose the probability density function (pdf) of \b\ is 
[x) and the cumulative distribution function (cdf) is -Fj^| (x). Then the pdf of the k-ih. smallest \ h\ out of n 
independently and identically distributed (i.i.d.) |&|'s is 

As K oo, f^^{x) converges to a simplified asymptotic expression according to the following lemma 
EQl. 

Lemma IV.2: Let f'^^\x) be the pdf of the fc-th smallest element out of n i.i.d. random variables whose 
pdf's are f{x) and cdf's are F{x). As n — > oo, f^''\x) S{x — where S{.) is the Dirac delta 

function. 
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With Lemma rV.2, we can prove the following theorem, which provides sufficient and necessary condition 

for unreliable decision feedback in asymptotic sense. 

Theorem IV.3: As ^ oo, (n, G R if and only if \bn{i)\ < Fr^(«;), almost surely. 
Proof: When ^x< F7\k), 

■iKM 



P{{n,i) e R) = ^ P (\bnii)\ is the p-th smallest] |6„(i)| = 



p=i 



kKM p (^hn{i) \ = x\\hn{i)\ is thep-th smallest^ P {\hn{i)\ is the p-th smallest^ 



p 



(\k{i)\=x) 



^6{x-F^ 
f\b\{x) 

1. 



dr. as K ^ CO 



When (n, i) e R, 

/ _ X E^f r P(x< (k) 1 1 6„ (i) I is the p-th smallest) P f 1 6„ (i) | is the p-th smallest) 
p(.>F-(.)|(n,.).R) = ^ PiM.K)^ ^ 

Epfr/^-i(«)/|£|(a;)c;a; 
Piin.i) e R)KM 



= 0. 



This completes the proof. ■ 
With Theorem IV.3 we can obtain the asymptotic pdf f^{x) of the unreliable set constructed by the 



confidence metric. 



where u{x) is the unit step function, k = F-Hk) and C = g ^ is a normalizing constant. 

I I Jo J\b\(^) 

However, what we are really interested in is the distribution of the decision feedback error in the unreliable 
set. Let < x < 2, then 



P(|A6|<a;|R) = P (^\Ab\ < x \ \b\ < nj 
P {\/\h\ < X, \b\ < Kj 
\b\ < k) 
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{P {0 < Ab < x,\l - Ab\ < k) + P {-X < Ab < 0, \ - 1 - Ab\ < k)) 



(P(l - K < Ab < min(x, 1 + k)) + P(max(-x, -1 - k) < Ab < -1 + k)) 



^P {1 - k < \Ab\ < x) , ifl-k<x<l + k 



< 1 



if x > 1 + K, 







if X < 1 — K 



V 



Thus the asymptotic pdf of the decision error in R is 



/|Afe|(^) = -:f\Ab\{x)u{x - 1 + k)u{-X + K + 1). 



(38) 



With the same manipulation, we can obtain the asymptotic pdf of the errors in the reliable set R: 




1 



f\Ab\{x) {u{x - I- k) + u{-x - k + 1)) . 



1 - K 



(39) 



V. Asymptotic Performance Analysis 



A. Asymptotic performance analysis 

To evaluate the system performance of RS EQ-MUD, e.g. the bit error rate, Monte Carlo simulations 
can be used at the cost of a large amount of computation, especially when either the number of users or the 
search width is large. As an alternative, asymptotic analysis in the large system limit {K — > oo, — * oo 
while keeping ^ = (3) can be applied in view of the recent success of such analysis in evaluating individual 
optimal MAP multiuser detection (lO-MUD) lEIll 0. 

For simplicity, we consider only the case of L = Li = 2, as introduced in Section IV.E, and focus 
on the detection of the t-th symbol of user k. Since the power of the desired signal is distributed in two 
successive symbol periods, we can consider symbol bk{t) as being transmitted through two independent 
virtual channels, detected independently and combined to obtain the soft input to the decoders. In each 
channel, there exist two groups of interferers, each of which contains nK users: one from the MAI in the 
current symbol period and the other from the ISI in the previous symbol period. After despreading, the 
powers of the desired signal, interference and noise, which are denoted by {Qi]^^^ p ^=0 i 

{N{i)}^^^ ^, respectively, are illustrated in Fig. 3, where channel and 1 denote the symbol period t and 
t + 1, respectively. We can treat the ISI as MAI from kK virtual users, thereby unifying the ISI into the 
framework of MUD. 



16 



Due to the assumption of exponential decaying variance of channel gains in (8), we can obtain the param- 
eters in the equivalent channel with explicit expressions, which are given by 



E{\pti^,m 



-A(.+j)(i_e~2A)(gA_;^)2 
2A 



-^^■(e^- 1) 



1 - 



A 2A 



A 



2A 



if i > and j > 
if 2 = and j > 
if z = j = 



(40) 



E{\pl^ih^)\'} 



-^(i-i)(l_e"^)2 



if z = 



if z > 



(41) 



and 



Nit) = al + il- K)pE{Ab^} iC{i, i) + C(0, 1)) 



0,1, 



(42) 



where P = ^ is the system load. 



It is shown in ETl that in the large system limit, the multiaccess channel i, i = 1, 2, is equivalent to 
a single user AWGN channel with input signal power r]iQ{i) and noise power N{i), where rji is the (non- 
asymptotic) multiuser efficiency of the virtual channel i. Using the replica method developed in the statistical 
mechanics of spin glasses lfT4l . the multiuser efficiency rji can be obtained by solving the following equation 



0: 



1 

Vi 



j=0 ' '^'^ ^ 



(43) 



^^^^ and the expectation is taken over the distribution of soft 



tanh^ [z,/YijVi + lijVi) 

where Dz = ^ exp (^-^j dz, -fij = §g 
decision feedback, which can be obtained by Theorem IV. 3. 

On obtaining the multiuser efficiencies of both virtual channels, we can compute the multiuser efficiency 
after combining the results from both channels by 

y^ l ViQii) 



V 



A=0 N{i) 

T 



(44) 



By obtaining the distribution of the soft decoder output by simulation, the evolution of r] with respect to 
the iteration stages can be achieved. This results in the evolution of bit error rate as a byproduct. 



17 

B. Parallel interference cancellation 

If we set /v = 0, the RS EQ-MUD degenerates to PIC, which treats all decision feedbacks as reliable and 
eliminates the time consuming trellis search at the cost of some performance degradation. The multiuser 
efficiency of PIC is given by 

-spi Q{ii 

r? = fl!ziim. (45) 
VI. Numerical Results 

A. Bit error rate 

Figure 4 shows the bit error rate for different values of j^, where E^, denotes the energy per information 
bit and A^o — 2(T^ is the noise one-sided spectral density, under the conditions K — k — 0.1 and 13 — 1. 
The iteration times are labelled near the corresponding curves. The channel codes are assumed to all be the 
convolutional code (23, 33, 37)8 with constraint length 5 here, and in subsequent simulations. 

We can see that the bit error rate diverges when ^ is less than 2.5dB. If ^ is larger than this threshold, 
each iteration improves the performance, which converges to the single user performance with moderate 

B. Validity of the asymptotic analysis 

Figure 5 shows the comparison between the asymptotic analysis and simulation results with different 
system loads, where ^ = 4dB and the other configurations are the same as in Fig. 4. We can see that the 
asymptotic analysis matches the simulation results of the finite case quite well, even when the unreliable set 
is small (here kK — 3). Thus in the future experiments, we apply the asymptotic analysis for computational 
efficiency. 

C. Performance with different search widths and system loads 

Figure 6 shows the bit error rate with the system load (3 ranging from 0.8 to 1.2 and with the search 
width K — Q (namely, PIC), 0.2 and 1 (namely, optimal MAP EQ-MUD). The iteration times required for 
convergence are labelled on the figure. We can see that PIC achieves approximately the same performance 
of RS EQ-MUD with k = 0.2, at the cost of more iterations. Both PIC and RS EQ-MUD with k = 0.2 
diverge when P > 0.95, whereas optimal MAP EQ-MUD attains the single-user performance for almost all 



18 

system loads with slightly increasing iteration times. This means that the use of small search widths may 
incur considerable performance loss when compared with the optimal MAP EQ-MUD. 

Figure 7 shows the bit error rate with various search widths n in the two cases of/?= 1.1, 1^ = 4dB and 
j3 — 1.2, jj^ — 5dB. The required iteration times are labelled on the figure. We can observe the waterfall 
phenomenon near k = 0.5 and the iteration times are reduced as the search width increases, at the cost of 
rapidly ascending computation. The evolution of multiuser efficiency 77 with P — 1.1 and ^ = AdB, is 
given in Fig. 8, where 77 increases to more than 0.9 with different required iteration stages when k > 0.5. 

VII. Conclusions 

Mitigation of ISI and MAI is of considerable importance to the performance of DS-CDMA systems. We 
have formulated the optimal MAP EQ-MUD for such systems with channel codes. To reduce the prohibitive 
computational cost of MAP EQ-MUD, we have proposed the RS EQ-MUD by decomposing the trellis 
and classifying the decision feedback into reliable and unreliable sets with confidence metrics. Asymptotic 
analysis has been used to evaluate the system performance and is shown to match simulation results for 
finite numbers of users. Numerical results show that the RS EQ-MUD can achieve close to single-user 
performance with moderate With a reasonable search width, RS EQ-MUD outperforms the conventional 
PIC with higher user capacity and fewer required iterations. 
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Illustration of the two equivalent independent channels for performance analysis 
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Bit error rate versus system load 
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Fig. 7 

Bit error rate versus search width 
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Multiuser efficiency at different iteration stages 



